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Abstract 

Defect-induced ferromagnetism has triggered a lot of investigations and controversies. The major 
issue is that the induced ferromagnetic signal is so weak that it can sufficiently be accounted for 
by trace contamination. To resolve this issue, we studied the variation of the magnetic properties 
of SiC after neutron irradiation with fluence covering four orders of magnitude. A large param¬ 
agnetic component has been induced and scales up with defect concentration, which can be well 
accounted for by uncoupled divacancies. However, the ferromagnetic contribution is still weak and 
only appears in the low fluence range of neutrons or after annealing treatments. First-principles 
calculations hint towards a mutually exclusive role of the concentration of defects: a higher defect- 
concentration favors a larger magnetic interaction at the expense of spin polarization. Combining 
both experimental and first-principles calculation results, the defect-induced ferromagnetism can 
probably be understood as a local effect which cannot be scaled up with the volume. 


"Electronic address: s.zhou@hzdr.de 
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I. INTRODUCTION 


The magic of magnetism was disclosed in the early 20th century with the development 
of quantum mechanics. The Heisenberg model has since then been extremely successful to 


understand magnetism and magnetic materials |l|. All previously identified ferromagnetic 
bulk materials contain elements with partially filled 3d (4d) or 4/ (5/) shells. A fundamen¬ 
tal question is whether materials containing only s or p electrons can be ferromagnetic. For 
nearly two decades, there have been various theoretical and experimental studies devoted 
to clarifying this question j2-2Q|. It turns out that materials with completely hlled 3d or 4/ 
shells or with only s or p electrons can be ferromagnetic when they contain defects. Among 
those materials, graphite/graphene and oxides attract particular attention due to experi¬ 


mental evidence reported by various groups 


t 5, 10, 21 


27l |. However, the experimentally 


measured ferromagnetism remains a weak signal slightly above the detection limit of sen¬ 


sitive SQUID magnetometry 


jg, 17, 28 


3l| . The very weak magnetization not only limits 


the practical applications, but also raises questions on the fundamentals of defect induced 
ferromagnetism. On the one hand, measurement artifacts in SQUID magnetometry may 


occur: improper mounting o 


ferromagnetic like signal 


32 


samples and wrong use of sample holders can easily generate 


34| . On the other hand, the debate over the purity of graphite 


and oxide substrates continues in parallel 


35 


39| . Pristine graphite and oxides substrates 


are often ferromagnetic due to different contaminations or due to intrinsic defects 


36 


40 


42|, which hamper the interpretation of the observed ferromagnetism. Indeed, Nair et al. 


reported that there is only paramagnetism in graphene after introducing adatoms or defects 


43|. 


Therefore, the understanding of defect-induced ferromagnetism in materials without par¬ 
tially hlled 3d or 4/ electron shells is far from satisfactory. It is rather calling for an 
investigation in which the following requirements should be fulhlled: 


• The pristine materials should be well controlled with the highest purity grade possible. 

• The materials should be supplied in a large quantity with identical properties, such 
that one can perform a series of experiments using identical specimens. 

• The materials should be free of elements with 3d or 4/ electrons. 
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• The induced effect should exist in a large volume to measure a large enough magnetic 
signal. 

The last point was firstly proposed by Coey et ah, who suggested to measure a bulk 
g.aphite-„odule to goto a large fertontagnetie etgnal g However, arator antounte of ntag- 
netite, kamacite, etc, also appear in their graphite samples [2[. These secondary phases are 
responsible for about two-thirds of the observed magnetization and the remaining one-thirds 
is attributed to the graphite-nodule. 

After screening by these facts, Si and SiC are the best candidates. Defect-induced ferro¬ 
magnetism was revealed in SiC after neutron and ion irradiation 


30 


31| . Recently, it was 


suggested that the p electrons from the carbon atoms are mainly responsible for the long- 
range ferromagnetic coupling in SiC after ion irradiation which is similar to the obser¬ 


vation in graphite 


45| . In this work, we start with 4H-SiC, high purity semi-insulating SiC. 


Neutron irradiation was used to introduce defects in the whole sample volume over a large 
range of defect concentrations. The large “magnetic volume” and well controlled pristine 
materials allow for reliable interpretation for the following experimental observations: (1) 
paramagnetism increases with neutron fluence and saturates at the end; (2) ferromagnetism 
only occurs in a narrow fluence range although it is weak; (3) thermal annealing param¬ 
agnetic defective SiC can induce weak ferromagnetism. Together with density-functional 
theory (DFT), we show that there is an intrinsic limit to increase the ferromagnetic contri¬ 
bution to the magnetization. 


II. EXPERIMENT 

Commercial semi-insulating 4H-SiC single crystal wafers (Cree) were used for our experi¬ 
ment. Particle induced X-ray emission (PIXE) was performed to check magnetic impurities 
in the pristine sample. The concentrations of magnetic impurities (Fe, Co and Ni) are 
proved to be below the detection limits (shown later). Neutron irradiation was performed at 
chambers DBVK and DBVR at the BER II reactor at Helmholtz-Zentrum Berlin at a tem¬ 
perature less than 50 °C. The fluxes of DBVK are 1.08x10^^ cm“^s“^ of thermal neutrons, 
7.04x10^^ cm“^s“^ of epithermal neutrons, and 5.80x10^^ cm“^s“^ of fast neutrons, respec¬ 
tively and those of DBVR are 1.50x10^^ cm“^s“^, 6.80x10^^ cm“^s“^, 5.10x10^^ cm“^s“\ 
correspondingly. Comparing with epithermal or fast neutrons, thermal neutrons only pro- 
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herefore, only epithermal and fast neutrons are accounted 
46| . The neutron fluence spanned a large range from 3.28x10^® 


duce negligible displacement.^ 
for in the fluence calculation 
cm“^ to 3.50x10^® cm“^ covering all possibilities from light damage to near amorphization. 


It is worthy to note that the applied neutron fluence covers the range used in Ref. 


^1 


and is extended to much higher fluence values by two orders of magnitude. For structural 
characterization, /i-Raman spectroscopy has been performed by using a Nd:YAG Laser with 
532 nm wavelength in the scattering geometry using a liquid nitrogen cooled charge-coupled 
device camera. Positron annihilation Doppler broadening spectroscopy (DBS) was applied 
to clarify the nature of defects using the mono-energetic slow position beam ’’SPONSOR” 


47l |. High-resolution transmission electron microscopy (HR-TEM) investigations were done 


with an image-corrected Titan 80-300 microscope (FEI) operated at an accelerating voltage 
of 300 keV. The magnetic properties are measured by a superconducting quantum interfer¬ 
ence device magnetic property measurement system and vibrating sample magnetometers 
(SQUID-MPMS and SOUID-VSM, Quantum Design). The electron spin resonance (ESR) 
spectroscopy was performed at 9.4 GHz using a Bruker spectrometer (Bruker ELEXSYS 
E500). 

III. RESULTS 


A. Virgin vs. irradiated SiC 

Before we start to describe the detailed results, we hrst present the key facts for the virgin 
SiG substrates employed in our investigation. 

To check the trace presence of transition metals (Fe, Go and Ni, etc.) in our SiG sub¬ 
strates, we have performed PIXE using 2 MeV protons with a broad beam of around 1 mnY, 
as shown in Figure [H PIXE is a sensitive method to detect trace impurities in bulk volume 
without signihcant structural destruction [35|, l39| . In the spectrum, the sharp peak is the Si 


K-line X-ray emission. The broad bump is due to the secondary electron Bremsstrahlung 
background. For commercially available, purest graphite, there is detectable transition metal 
contamination (mainly Ti, V, Fe, Ni) as revealed by PIXE 39|. In sharp contrast, the tran¬ 
sition metal impurity in the SiG we used for this study, if there is, is below the detection 
limit of around 1 ppm (Fig. [T]). 
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FIG. 1: PIXE spectrum for the 4H-SiC wafer by a broad proton beam. Within the detection limit, 
no Fe, Co or Ni contamination is observed. 

Figure [2] shows a comparison of the magnetic properties measured at 5 K between the 
virgin SiC and a specimen irradiated with neutrons with a relative low fluence value of 
4.68x10^^ cm“^. The results are presented without any background subtraction but only 
normalized to the mass of the samples. The difference between the virgin and irradiated 
SiC is obviously not trivial: the virgin SiC shows purely diamagnetic behavior, while the 
irradiated one shows a paramagnetic component as revealed by the large deviation from the 
linear dependence on the magnetic held. The inset of Fig. [2] shows the zoom of the low 
held part. Besides the paramagnetic component induced by neutron irradiation, a small 
ferromagnetic contribution also shows up with its saturation magnetization around 1% of 
the paramagnetic signal. The defect-induced paramagnetism and ferromagnetism are exactly 
the key facts we are going to discuss in this manuscript. 


6 










Field (Oe) 


FIG. 2: Comparison of the magnetic properties: virgin vs. neutron irradiated SiC. Inset: the zoom 
of the low field part. Neutron irradiation induces signihcant, unambiguous magnetic variation in 
SiC. 

B. Structural properties 


1. Raman spectroscopy 


The damage to the crystallinity of SiC upon neutron irradiation is verihed by Raman 
spectroscopy. Figure [3] exemplarily shows the Raman spectra for the virgin sample and 
the samples irradiated with different neutron fluence values as indicated in the hgure. The 
typical Raman scattering modes for 4H-SiC are resolved: folded transverse optic (FTO) and 
longitudinal optic (FLO) modes j^|. Note the different scale in the y-axis: with increasing 
fluence, the intensity of the Raman scattering modes is dramatically reduced, as observed 
in Refs. [30, 


31| . This behaviour directly reflects the increasing disorder of the crystalline 


material. For the sample with largest fluence, the Raman modes are much broad and hardly 
resolved due to the large amount of defects l49|. However, after annealing at 900 °C for 
15 min the crystalline order can be substantially recovered. 
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FIG. 3: Raman spectra for virgin and neutron irradiated 4H-SiC single crystals. The folded 
transverse optic (FTO) and longitudinal optic (FLO) modes are identified. With increasing neutron 
fluence the Raman scattering intensity is decreased and the peak is broadened. 

Figure m shows the Raman spectra of an irradiated SiC (fluence: 3.50x10^® cm“^) after 
annealing at different temperature. The annealing was performed in N 2 atmosphere at 800, 
900 and 1000 °C for 15 min. With increasing annealing temperature, the peaks of the SiC 
Raman modes become sharper and their intensity gradually increases. The increase of the 
peak intensity is related to the recovery of the crystalline lattice with annealing. However, 
a complete recovery of the lattice is not achieved. It has been shown that after annealing at 
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FIG. 4: Raman spectra for neutron irradiated and thermally annealed 4H-SiC single crystals. The 
annealing was performed in N 2 atmosphere 15 min at the temperature indicated. With increasing 
annealing temperature, the structural damage has been healed. 

1450 °C for 10 min the Raman peaks are comparable to the pristine SiC 



2. Positron annihilation spectroscopy 


Positron annihilation Doppler broadening spectroscopy (DBS) is an excellent technique 
to detect ^en volume defects from clusters consisting of several vacancies down to a mono¬ 


vacancy 


47l |. The positron in a crystal lattice is strongly subjected to repulsion from the 


positive atom core. Because of the locally reduced atomic density inside the open volume 
defects, positrons have a high probability to be trapped and to annihilate with electrons 
in these defects by the emission of two 511 keV photons. The Doppler broadening of the 
511 keV annihilation line is mainly caused by the momentum of the electron due to the 
very low momentum of the thermalized positron. The Doppler broadening parameter S, 
obtained from the 511 keV annihilation line, reflects the fraction of positrons, annihilating 
with electrons of low momentum (valance electrons). In this study, the S parameter is 
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defined as the ratio of the counts from the central part of the annihilation peak (here 510.17 
keV - 511.83 keV) to the total number of counts in the whole peak (498 keV - 524 keV). 
Therefore, the S parameter is a measure for the open volume in the material. It increases 
with increasing size of the particular open volume defects. 

According to our previous positron annihilation experiments 


30 


31| . divacancies have 


been identihed in neutron or ion irradiated SiC. Figure [5](a) shows the measured S parame¬ 
ters versus the incident positron energy. The plateau of the S parameters of the irradiated 
samples above a positron energy of 2 keV corresponds to positron annihilation in the defects 
created by neutron irradiation. Defects are homogenously distributed along the depth as 
expected for neutron irradiation. The S parameters of irradiated samples are larger than 
that of the virgin SiC. S does not increase signihcantly with increasing neutron fluence, 
but does increase after annealing. Figure |5](b) displays the S — W plot for different sam¬ 
ples. A straight line is obtained in the S — W representation for different SiC samples. It 
shows that the same type of defects, i.e., open volume damage {m, 3; exist in different SiC 
samples. A relation between the ratio Sdefect/Sbuik and the number of agglomerated VsjVc 


divacancies in 6H-SiC is published as a scaling curve in Ref. ImI As a rough estimation, the 


post-irradiation annealing at 900 °C led to a defect agglomeration up to a size of 8 Vsi^c 
divacancies. 

One has to note that PAS is only sensitive to negatively charged or neutral open volume 
defects. The S parameter increases with the size of the open volume defects, but not 
necessarily with the amount of defects. Raman and magnetization are measurements of 
total amount of defects. This explains why the S parameter increases after annealing, 
indicating coalescence of voids to larger but fewer ones, but the Raman and magnetization 
measurements show the decrease of amount of defects. 


3. High-resolution transmission electron microscopy 

With the aim to visualize the defects, we performed high resolution transmission electron 
microscopy on selected samples: a pristine 4H-SiC, 431-56 (irradiated with neutron at the 
fluence of 3.50x10^® cm“^) and 431-56 after annealing at 900 °C for 15 min. Figure [6] dis¬ 
plays corresponding cross-sectional high-resolution transmission electron micrographs taken 
in [100] zone axis geometry. It is observed that the 4H-SiC stacking sequence (abcb) perpen- 
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FIG. 5: (a) Mean S parameter vs. incident positron energy for virgin, neutron irradiated (431- 
50 with fluence of 3.50x10^® cm“^ and 431-56 with fluence of 3.50x10^® cm“^) and thermally 
annealed (431-56) 4H-SiC single crystals. The annealing was performed in N 2 atmosphere at 900 
°C for 15 min. (b) S' — IT plot of DBS results measured on different SiC crystals shown in (a). 
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- See more at the original publication Phys. 


Rev. B 92, 174409 (2015) due to the size limit. 

FIG. 6: Cross-sectional HR-TEM micrographs of 4H-SiC single crystals: (a) Pristine SiC (b) 431- 
56 (irradiated with neutron fluence of 3.50x10^® cm“^) and (c) 431-56 after annealing at T = 900 
°C. 


dicular to the surface has been kept well during neutron irradiation up to the applied fluence. 
There are no observable contrast differences between the samples, even after annealing. This 
hnding is due to the fact that neutron irradiation (in relatively low fluence regime compared 
with literature) mainly creates point defects. According to the magnetization results shown 
later, the average defect concentration for the sample with the largest fluence is below 0.1%. 
Furthermore, it should be kept in mind, that classical TEM specimen preparation by sawing, 
grinding, polishing, dimpling, and final Ar ion milling may also lead to the generation of 
point defects, resulting in minor contrast variations in the HR-TEM micrographs compared 
to a completely defect-free single crystalline 4H-SiC specimen. Therefore, TEM seems to be 
not the right tec hnique t o clearly visualize the defects in the relatively low fluence regime. 

)|, defects appear to be visible by TEM when the neutron fluence 


As shown in refs. 


491, 


53-55 


is above 10^^ cm i.e. 100 times larger than in our case. 
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FIG. 7: (a) Magnetization of all irradiated samples measured at 5 K as a function of field [(M(H)\, 
(b) Fits of the magnetization measured at 5 K for the sample with a fluence of 3.50x10^® cm“^ 
using the Brillouin function with different values of J. 

C. Paramagnetism 

Figure [71(a) presents the magnetization [(M(HJ\ of neutron irradiated SiC after subtract¬ 
ing the diamagnetic background which is measured for the pristine sample at 300 K. It is 
clearly seen that all samples exhibit a paramagnetic component. The magnetization in¬ 
creases with increasing neutron fluence. The paramagnetism can be accurately described by 
Brillouin function: 


13 










71 ^/ \ 717- T r2t/ -|- 1 ,2(7 -|- 1 . 1 1 / ^ \1 

M{a) = NJfiBg[^j—coth{^j—a) - ^coth{—)] 


( 1 ) 


where the g factor is about 2 obtained from electron spin resonance measurement (shown 
later), fiB is the Bohr magneton, a = gJfiBH/kBT, kB is Boltzmann constant and N is the 
density of spins. As exemplarily shown in Fig. [71(b), the data can be htted nicely by J = 
1. Moreover, a value of J larger or smaller than 1 results in a pronounced deviation from 
the experimental data. This htting allows us to conclude that the observed paramagnetism 
is due to a single species with moment jj, = gJjJ-B = 2/1^. This conclusion is in excellent 
agreement with our first-principles calculation (shown later): the induced magnetism is due 
to VsjVc divacancy carrying a magnetic moment of 2 /i^. As a comparison, the induced 
magnetism in graphite is generally attributed to single vacancy 5fil-l59| . In ZnO, the origin 


of the local moment is more complicated. Bot 
be responsible for the induced magnetism 


1 Zn and O vacancies have been suggested to 


60 


62l |. On the other hand, oxygen atoms on 


the polar ZnO surface can be spin polarized and undergo a long-range magnetic interaction 
63|-|65|. 

Figure IH] shows the temperature dependent magnetization [M(TJ\ for samples with dif¬ 
ferent neutron fluences under an applied held of 10000 Oe. As expected, one contribution 
represents the diamagnetic background, which is essentially temperature independent and 
dominates at higher temperature. The Curie-like paramagnetic component shows a strong 
temperature and neutron huence dependence. In a self-consistence, the M(T) curves can be 
well htted according to the Curie law (Eq. [2]) by using the same J and N obtained from 
the corresponding M(H) curve htting: J = 1 and N = 1.66xl0^®/g for the sample shown 
in Fig. [71(b) and Fig. [8](b). 


M _ J{J+l){gi,Br ,,1 

H ikaT ' ’ 

Figure [9] depicts the magnetization of irradiated SiC after thermal annealing at diherent 
temperature. As shown in the Raman investigation (Fig. [1|), thermal annealing heals par¬ 
tially the crystalline damage. After annealing, although paramagnetism is still the major 
component in the magnetic signal, the magnetization is much reduced. The removal of point 
defects results in the reduction of paramagnetism. The appearance of the ferromagnetic cou¬ 
pling after annealing will be shown in the next subsection. 
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FIG. 8: (a) The magnetic moment of all irradiated samples measured at 10000 Oe as a function of 
temperature \M(T)]. (b) For the sample irradiated with the largest fluence: the black symbols are 
experimental data and the red solid curve is the fitting result by equation (2). 

Paramagnetic centers in the neutron irradiated SiC have also been detected by ESR. As 
shown in Figure [TOl in the broad field range, there is only one sharp resonance peak due 
to paramagnetic electrons. Within the detection limit, there is no broad resonance peak, 
which could be related with ferromagnetic resonance. The inset shows a narrow scan with 
smaller field step at the paramagnetic resonance peak. The ESR spectrum exhibits a central 


15 






















FIG. 9: Magnetization of irradiated SiC after annealing at different temperature: measurements 
at 5 K as a function of field [(M(H)]. 


line and some weak hyperfine lines. The g factor is calculated to be around 2.005, which 
is the characteristic of free electrons. From the line shape, the defect is very probably the 


divacancy in 4H-SiC 
et al. 


66| . The hyperhne structure is not well resolved as in the work by Son 


66| . which is probably due to the crystal degradation upon neutron irradiation. 


D. Ferromagnetism 


In Refs. 


30 


3l| . a sizeable ferromagnetism was observed for neutron or ion irradiated 


SiC. Note that in Ref. 


30| . the maximum neutron fluence (if counting only the fast neutrons) 


is around 5.60x10^^ cm“^. The authors show magnetization measurements at low held and 
a magnetic hysteresis loop was observed. We also checked our samples carefully whether 
a hysteresis appears in the low held range. Interestingly, besides the large paramagnetic 
component we also observe a ferromagnetic hysteresis in samples with huence of around 
5.60x10^^ cm“^. However, the hysteresis is not resolvable when the huence is higher than 
6-7x10^^ cm“^. In Fig. [TTl we show the detailed magnetization measurements for three 
representative samples: sample 431-50 with an intermediate huence 4.68x10^^ cm“^, and 
sample 431-56 with the highest huence of 3.50x10^® cm“^ before and after annealing. Sam- 
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FIG. 10: ESR spectrum for a 4H-SiC sample irradiated with neutrons with a fluence of 8.38x10^® 
cm“^. The field is applied perpendicular to the c-axis. The inset shows the detailed hyperfine 
structure (indicated by arrows) of the resonance peak. 

pie 431-50 shows a clear hysteresis added on the paramagnetic component (at 5 K) or on the 
diamagnetic backgronnd (300 K). In contrast, for sample 431-56, at both 5 and 300 K, there 
is no hysteresis resolvable and the paramagnetism dominates at low temperatnre. However, 
after annealing at 900 °C for 15 min, sample 431-56 shows a ferromagnetic component at 
5 and 300 K. The paramagnetic component is drastically rednced to below 10% compared 
with the non-annealed sample dne to the redaction of defects as confirmed by Raman spec¬ 
troscopy (Fig. |3]). Moreover, we also observed a similar ferromagnetic component in several 
samples with the fluence of the order of 2-6x10^^ cm“^ as well as in 6if-SiC samples with 
a large fluence of 3.50x10^® cm“^) after annealing. The saturation magnetization for the 
ferromagnetic component in different samples is in the range of 1-5x10“^ emu/g and is only 
around 1% of the paramagnetic component. 

The appearance or disappearance of the weak ferromagnetic component is confirmed 
by zero-field-cooled and field-cooled magnetization measurement (ZFC/FC) shown in Fig. 
[TTT d. e, f). The ZFC magnetization was measured by cooling the sample from 350 K to 5 
K with zero field, then a field of 100 Oe was applied and the magnetization was measured 
during warming up. The FC curves were measured by cooling the sample in a field of 100 
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Oe during cooling. This approach is often used to verify if the measured specimen is fer¬ 
romagnetic, paramagnetic or superparamagnetic. If the specimen contains a ferromagnetic 
component and if the held during measurement is smaller than its coercive held, the ZFC/FC 
magnetizations will show an irreversibility with the FC magnetization larger than the ZFC 
magnetization. If the specimen contains only paramagnetism, the ZFC/FC magnetization 
should overlap with each other. For samples that contain superparamagnetic components, 
a blocking behavior at low temperature, i.e. an increase in the ZFC magnetization with 
temperature, will appear. For samples 431-50 [Fig. [TTT dl] and 431-56 after annealing [Fig. 
mi)], a diherence is observed up to 300 K in the ZFC/FC magnetization, while for sample 
431-56 the FC curve superimposes the ZFC curve. Moreover, the data did not indicate 
any blocking temperatures that can be associated with superparamagnetic behavior, which 
was reported for ion implanted SiC [^. Figure [12] shows the magnetization vs. field at 
the low field range for samples after annealing at different temperature. Annealing at rel¬ 
atively low temperature (in contrast to 1450 °C normally used for SiC 50|) leads to the 
appearance of ferromagnetic hysteresis. However, the ferromagnetic component is clearly 
snppressed when the annealing temperatnre is increased to 1000 °C. As shown in Fig. |9] 
the paramagnetism in the sample annealed at 1000 °C is much reduced in comparison with 
samples annealed at 800 or 900 °C. Apparently, annealing at temperatures higher than 900 
°C leads to the annihilation of both isolated and agglomerated divacancies. Moreover, a 
similar annealing behavior (the reduction of paramagnetism and the appearance of a weak 
ferromagnetic component) has been observed in 4H-SiC with lower fluences as well as in 
6H-SiC (not shown). 

As shown by our comprehensive investigation on the magnetic properties of 4H-SiC npon 
neutron irradiation with a much broader fluence range than that used in Ref. (espe¬ 
cially in the large fluence range), the paramagnetic component scales up with fluence (the 
concentration of defects), but the weak ferromagnetic contribution appears only for samples 
with relatively low fluence or after annealing. Figure fldf a) shows the fitted paramagnetism 
depending on neutron fluence. Those samples denoted by the red circles and by the blue 


triangle appear to have a ferromagnetic component. As shown in Refs. 30|, 1^, the divacan¬ 


cies in SiC can explain those magnetic properties. The scaling between the concentration 
of paramagnetic centers and the neutron fluence provides solid evidence that the defects 
created by irradiation are directly responsible for the local moments. However, the long- 
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FIG. 11: Magnetization vs. field at the low field range for samples with fluences (a) 4.68x10^^ 
cm“^, (b) 3.50x10^® cm“^ and (c) 3.50x10^® cm“^ after annealing and ZFC/FC magnetization 
for the same set of samples (d-f). 


range-coupling between the local moments is not always favorable. The ferromagnetism 
only appears when the defect concentration (represented by the density of paramagnetic 
center and by the neutron fluence) is in a narrow window (as indicated by the grey bars 
in Fig. ITSf ai.i This is consistent with the ferromagnetism observed in ion implanted SiC 
0, 68|. Nevertheless, the saturation magnetization for the ferromagnetic component in dif¬ 
ferent samples is much smaller than the paramagnetic component. In the next section, we 
employ hrst-principles calculations to understand why the defect-induced ferromagnetism is 
weak. 


E. Why is the ferromagnetic signal weak? 


According to our previous positron annihilation experiments 30|, 1^, divacancies have 


been identihed in neutron or ion irradiated SiC and are the most likely origin for the mea¬ 
sured ferromagnetism. The ferromagnetism only occurs in samples with relatively low neu- 
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FIG. 12: Magnetization vs. field at the low field range for samples after annealing at different 
temperature. Neither the diamagnetic background nor the paramagnetic component has been 
subtracted from the measured signal. 

TABLE I: Magnetic coupling between V^iVc-induced local moments in SiC. The coupling is anti¬ 
ferromagnetic between neutral or positively charged Ysi^c while ferromagnetic between negatively 
charged ones. As the arrangement direction of Vsi^c is along the c axis, the coupling almost does 
not exist. 


dv (A) c- 

or a-axis 

charge state Eafm-^fm ( 

meV) Jo (meV) 

10.11 

c axis 

(0, 0) 

-0.55 

-0.14 

12.33 

a axis 

(0, 0) 

-39.45 

-4.93 

12.33 

a axis 

(+1 +1) 

-251.35 

-125.68 

12.33 

a axis 

(-1, -1) 

98.73 

49.36 

15.42 

a axis 

(0, 0) 

-5.77 

-0.72 

15.42 

a axis 

(-1-1, 4-1) 

-179.81 

-89.90 

15.42 

a axis 

(-1, -1) 

33.31 

16.66 

18.50 

a axis 

(0, 0) 

-0.93 

-0.12 

18.50 

a axis 

(+1, -1-1) 

-49.75 

-24.88 

18.50 

a axis 

(-1, -1) 

7.43 

3.71 
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FIG. 13: (a) The fitted N (the density of paramagnetic centers) for samples with different neutron 
fluence {(j)) and the sample with the highest fluence after annealing at 900 °C for 15 min. The dashed 
line is only for guiding eyes to indicate the trend. The grey bars indicate the range of samples with 
a resolvable ferromagnetic component, regarding to N and neutron fluence, respectively, (b) N/cj) 
vs. neutron fluence. 21 













FIG. 14: (a) Spin-polarization energy AE^p as a function of the distance dy between the adjacent 
vacancies. AE^p increases with increasing dy and the critical value is around 10 A. (b) The 
ferromagnetic (or antiferromagnetic) exchange coupling energy between charged (-1, -1) or (-1-1, 
-1-1) divacancies as a function of dy. The corresponding supercell structures with one Vsi'^C are 
marked with the numbers and multiply signs. 


tron fluence or after annealing and is mnch weaker compared with the induced paramag¬ 
netism. These facts motivated us to perform some theory work to understand the puzzle. 
First-principles calculations were performed by using the Cambridge Serial Total Energy 
Package |69|. Spin-polarized electronic structure calculations were carried out using the 


Perdew-Burke-Ernzerhof functional 


70| for the exchange-correlation potential based on the 


generalized gradient a ppr oximation. The core-valence interaction was described by ultra- 
soft pseudopotentials [7l|. The cutoff energy for the plane-wave basis is set to 310 eV. 
Full optimization of the atomic positions and lattice parameters was carried out with con- 
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vergence threshold of the remanent Hellmann-Feynman force 0.01 eV/A. 4H-SiC supercells 
with various sizes containing one axial divacancy (V 5 jVc) of built to obtain the variation 
of spin polarization and the magnetic coupling depending on the distance between defects. 
Actually, the calculation results show no significant difference between 4H- and 6H-SiC. 

As a hexagonal structure, 4H-SiC has a lattice constant c (out-of-plane) different from a 
and b (in-plane), so the distance between adjacent Vsi^c divacancies along different axes 
is also different from each other. Usually, the magnetic moments couple with each other 
preferentially along the shortest distance, which will be chosen as the coupling distance. It 
is found that the coupling is along the a-h plane in a SiC supercell, which will be explained 
hereinafter. Therefore, the distance between adjacent V^iVc in the a-b plane is chosen as 
the coupling distance {dy) between them. The spin-polarizing energies AEgp (the energy 
difference between the spin-unpolarized and spin-polarized states) of supercells as a function 
of dy are shown in Fig. [141 It is found that AE^p increases when dy increases, indicating 
that the spin polarization is stable under low Vsi^c concentrations. When dy is larger than 
around 10 A, AE^p changes from negative to positive and the spin polarization becomes 
energetically favored. Here, the supercells with dy less than around 10 A are designed for 
obtaining the critical value of dy switching the sign of AE^p. Besides, each V^iVc yields 
a local moment of 2.0 when dy is more than 12.34 A. This result is in good agreement 
with our experimental data shown in Figs. [7] and [HI The experimental results show that 
the induced paramagnetism is due to a single species with a moment of 2 Moreover, as 
shown in Fig. [THf b). the normalized density of paramagnetic center by the neutron fluence 
(i.e. the measure of the amount of defects), N/cf) decreases with the neutron fluence. This 
trend again indicates that a low defect-density favors spin-polarization as shown in Fig. 
[UKa). 

In order to investigate the magnetic coupling (ferromagnetism or antiferromagnetism) 
between these V^iVc-induced local moments, we put two 4H-SiC 4x4x1 supercells side 
by side along the a- or along the c- direction. Each 4x4x1 supercell contains one VsNc- 
Note that these antiferromagnetic structures are designed only for obtaining the magnetic 
interaction. The energy difference between the antiferromagnetic and ferromagnetic phases 
is Eafm — EpM = 8Jo{dy)S'^ or Eafm — EpM = 4Jo(d\/)>S'^ for size doubling along a or c 
directions, respectively, according to the nearest-neighbor Heisenberg model, where Jo{dy) 
is the nearest-neighbor exchange interaction as a function of dy and S is the net spin of the 
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VsiVc states. In previous work, it has been shown that the intrinsic defects in high purity 


SiC can be neutral, positively or negatively charged 

n 


72| . The possible impurities, such as 
Boron and Nitrogen (see ref. 173|), with the concentration of the order of can 


also modify the charge state of V^iVc. Therefore, we calculate the interaction for different 
charge states. As shown in Table I, VsjVc does not couple with each other along the c 
axis, even when dy is as small as 10.11 A. The exchange interaction is much stronger along 
the a axis and strongly depends on the charge state. Neutron or positively charged V^iVc 
antiferromagnetically couples with a maximal exchange interaction of 125.68 meV when dy 
is 12.33 A. The result is also consistent with the previous report 30|, in which the coupling 
between neutral VsjVc is also antiferromagnetic. Negatively charged V^iVc divacancies 
favor a ferromagnetic coupling. However, the exchange interaction dramatically decreases 
from 49.36 meV to 3.71 meV then dy is increased from 12.33 to 18.50 A. 

The mutual dependence of spin-polarization and exchange interaction on the concentra¬ 
tion of defects is not only suggested for SiC, but also for graphene/graphite [74 l75| as well 


as for nitrides 


16 


76| and oxides j^. To develop the long-range ferromagnetic coupling, it 


is also crucial to keep the stacking order in graphite |78|. 

In Figure im b). we plot the ferromagnetic exchange interaction energy depending on dy. 
As shown in Fig. m a) and (b), the ferromagnetism induced by divacancies requires that 
divacancies are negatively charged with dy in the range from 12.33 A to 15.42 A. The required 
dy would correspond to a concentration of divacancies of the order of 3-6 xlO^'^ cm“^ (l/dy). 
This concentration is around 1 at.% of SiC and well above the maximum concentration 
for the local moments we obtained by neutron irradiation (see Fig. [H- However, after 
irradiation with the largest neutron fluence the density of paramagnetic centers (A) starts 
to saturate probably due to the appearance of local regions which contains a large density 
of defects and tends to be non-spinpolarized as shown in Fig. [Mlf ab Moreover, it is also 
crucial to keep the crystalline ordering for long-range ferromagnetic coupling TSj- In this 
sense, it is unrealistic to expect ferromagnetic coupling throughout the whole bulk sample. 
In principle, it is in agreement with our experimental observation. The experimentally 
observed weak ferromagnetism probably can be understood as a local effect: only some 
particular regions with nm or pm dimension can accommodate a high concentration of 
divacancies. These regions form ferromagnetic bubbles with very strong interaction (see 
Table I), leading to the high Curie temperature. Presumably, these regions could be at the 
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surface, an interface or a grain boundary. However, according to our calculation the spin- 
polarization will be unfavorable if dy is too small. Therefore, the narrow fluence window of 
ferromagnetism shown in Fig. [12] as well as in neon implanted SiC Slj] can be speculated: 
In the beginning of neutron irradiation, the defect concentration increases, and some local 
regions accommodating a larger concentration of divacancies appear to be ferromagnetic. 
When dv in the local regions reaches the critical value of about 10 A, the increase of the 
defect concentration will suppress the spin polarization. There may appear some new local 
regions with the concentration of divacancies reaching to the level to induce ferromagnetism, 
but the damage to the crystalline structure by irradiation will weaken the coupling. Also, 
the existence of minor and major carriers in SiC can lead anti-ferromagnetic coupling and 
cancel the ferromagnetic coupling to some degree. Therefore, the magnetic moments can 
not be increased by increasing the fluence of irradiation all the time, while, the amount of 
paramagnetic centers can be scaled up with neutron fluence and show saturation at very 
large fluences. The locally accumulation of ferromagnetic defects has also been suggested 
for defective ferromagnetic oxides [l4, 3, 

To clarify if the defect induced ferromagnetism is associated with some local regions, 
it rather calls for a sophisticated investigation using methods which are sensitive to the 
magnetic order and provide an enough spatial resolution at the same time. Spin polarized 


scanning tunneling microscopy 


57( 1 or muon spin relaxation 80|, 1 m 1 could provide us a 


further insight, which is beyond the scope of the present paper. 


IV. CONCLUSION 

To investigate defect-induced magnetism in SiC, we applied a broad neutron fluence cov¬ 
ering four orders of magnitude to generate defects in a controlled way. A huge paramagnetic 
contribution is observed and scaled up with neutron fluence, whereas ferromagnetism only 
appears in a certain, rather low fluence range or after annealing treatments. First-principles 
calculations hint towards a mutually exclusive role of defects: A low defect concentration 
favors spin-polarization, but leads to negligible magnetic interaction. Combining both ex¬ 
perimental and hrst-principles calculation results, we conclude that defect-induced ferro¬ 
magnetism cannot be scaled up with the volume. It should be rather understood as local 
effect. Therefore, our investigation answers why the measured magnetization is low for 
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defect-induced ferromagnetism in SiC. 
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